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Abstract 

We show, both experimentally and theoretically, that sympathetic cooling of 87 Rb atoms in 
the \F = 2,m_F = 2} state by evaporatively cooled atoms in the \F = l,mF = —1) state can 
be precisely controlled to produce dual or single condensate in either state. We also study the 
thermalization rate between two species. Our model renders a quantitative account of the observed 
role of the overlap between the two clouds and points out that sympathetic cooling becomes 
inefficient when the masses are very different. Our calculation also yields an analytical expression 
of the thermalization rate for a single species. 
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Evaporative cooling has proved to be an efficient route toward Bose-Einstein condensation 
(BEC). However, not all the species are eligible for evaporation. For instance, the collisions 
may not be sufficient or may be forbidden as for fermions at low temperature B. In 
the case of a rare species, one may also want to avoid the large loss of atoms inherent to 
evaporative cooling. Sympathetic cooling is potentially a very good solution in those cases. 
It consists in putting the target sample we want to cool in thermal contact with a buffer 
gas that can be cooled to the desired temperature. This cooling can be done either by 
conventional cryogenics || |J, or by evaporative cooling as demonstrated by Myatt et al. 
||. These authors have produced dual BEC by sympathetic cooling of a gas of 87 Rb atoms in 
the \F = 2, mp = 2) state (target gas, noted F = 2 hereafter) by the use of an evaporatively 
cooled gas of 87 Rb atoms in the \F = l,mp = —1) state (buffer gas noted F = 1). In 
their experiment, the number of target atoms, although almost constant during the initial 
phase of the cooling, decreases significantly near to the condensation. In this paper, we 
investigate a situation similar to the one of Myatt et al, but with the number of atoms 
in the target gas constant in the last cooling stage, all the way down to BEC. This allows 
us to study quantitatively how the production of dual BEC relies on a precise control of 
the initial conditions. We also study the effect of the interspecies thermalization rate on 
sympathetic cooling efficiency, by changing the overlap between the buffer and the target 
gas. Our experimental observations agrees with the results of a calculation of the thermal 
contact between the atomic species. This calculation not only shows the crucial role of the 
overlap, but also why sympathetic cooling is very efficient for equal masses of the target and 
buffer gases. In addition, it allows us to recover in an analytic way a result previously found 
numerically for the thermalization rate of a single species. 

Our experimental apparatus has been described in a previous publication ||. After a 
laser cooling sequence, we transfer the 87 Rb atoms from a dark SPOT into the magnetic 
trap. During the transfer, an optical pumping pulse allows us to transfer most of the atoms 
in F = 1 (buffer), while keeping an adjustable fraction in F = 2 (target). Our iron-core 
Ioffe-Pritchard trap operates with a high bias field Bq between 50 and 250 Gauss, and the 
magnetic field curvature C along the dipole direction is linked to -Bo (the ratio B Q /C is 
fixed by construction to 1cm 2 ). After adiabatic compression, the magnetic field gradient in 
the strong (quadrupole) directions is G = lkG/cm. To cool down F — 1, we use standard 
rf forced evaporation. We characterize both species simultaneously by absorption imaging 
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after turning off the magnetic trap. During this release, the atoms experience a magnetic 
kick due to a Stern-Gerlach effect which allows us to spatially separate atoms in F = 1 and 
atoms in F = 2. We probe on the F = 2 — > F 1 = 3 transition after repumping the atoms 
from F — 1 to F — 2. With usual analysis, we derive the temperature of the two clouds 
from the image. 

In a first experiment, we have studied sympathetic cooling at a bias of 56 G, with initially 
10 8 atoms in F = 1 at a temperature of 300 /iK, and various initial numbers of atoms in 
F = 2. We observe no loss of target atoms during the whole cooling process 0. Indeed, 
the atoms in F = 2 see a trapping potential twice as steep as for F = 1, and as in ||, very 
few target atoms are evaporated at the beginning of the evaporation sequence. Moreover, at 
our high bias field, the non linear Zeeman effect prevents rf coupling of the \F — 2, mp = 2) 
trapping state to any non trapping state at the end of the evaporation || || . When starting 
with a very small (< 3 x 10 4 ) number of atoms in F = 2, we get a F = 1 condensate with 
typically 10 6 atoms. If we continue the evaporation, another condensate appears in F = 2 
at a lower temperature. For 3 x 10 4 atoms in F = 2, both gas condense simultaneously at 
a temperature of 200 nK, as shown figure |l|.a. If we start with a larger number of target 
atoms (but smaller than 8 x 10 4 ), a condensate appears first in F = 2 at T > 200 nK, and 
then in F = 1 at T < 200 nK. For a yet larger number of target atoms, we can obtain a 
condensate in F = 2 and no condensate in F — 1. If the number of atoms in F = 2 is more 
than 2 x 10 5 , BEC is observed in neither hyperfine state. 

These experimental observations can be explained with a simple model based on an energy 
budget. We assume that both species are always thermalized (see last part of this paper), 
and that the number N 2 of target atoms is constant. The total energy of the Ni buffer 
and N% target atoms in harmonic traps is therefore E = 3(Ni + A^/cb^- The buffer gas is 
evaporatively cooled with an energy cutoff rj k^T. We denote dNi the elementary (negative) 
variation of the number of buffer atoms. The corresponding (negative) energy variation is 
dE = dNi(r] + l)ksT 0. After rethermalization at a lower temperature T + dT, the total 
energy becomes E + dE = 3(Ni + dNi + A^2)A^(T + dT). Keeping only the first order, we 
have dT /T = a dNi/(Ni + N 2 ) with a = (rj — 2)/3. Assuming that r] is constant, we get: 
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In this formula, T im and iV{ m are the initial values, N 2 is constant and has been neglected 
with respect to iV{ m in the expression of T mm (as the temperature has to be reduced by 
several orders of magnitude, the initial fraction of atoms in the target gas must be very 
small). The temperature T mm is reached when all the buffer atoms have been evaporated 
(Aq = 0). To understand the different regimes experimentally observed, we first study the 
evolution of the phase space densities T>\ of the buffer and T> 2 of the target gas, in harmonic 
traps of mean frequencies u)\ and u 2 , during the evaporation. Equation |l| allows us to express 
T>i and T> 2 as functions of Aq only. Figure § represents the evolution of T>\ and T> 2 during 
the evaporation (while Aq is decreasing). We observe that the phase space density of the 
target gas V 2 always increases to reach its maximum value 

r max_ 2 -!7 ( NiBi a ^2 V () 

2 " Nt~ l) V 1 k B T^ ) {2) 
at the very end of the evaporation fill; this is because is constant and the temperature 



always decreases. We also observe that the buffer gas phase space density T>\ first increases 
|TT| to a maximum D™ ax that we write 



vr^vrM '^r" . (3) 



,lu 2 J (3a) 3a 

and then decreases to zero. At the intersection of the two curves, the phase space densities 
T>\ and T> 2 have the same value 

v^ = vr x ( 1 + (^) 3 ) 3a ■ ( 4 ) 

We note that a modification of the various parameters only changes the relative position of 
the two curves, in particular the ordering of T>^ = \ p max and D™^. This is what determines 
the outcome of sympathetic cooling. 

We specifically focus on the situation where T)^ is reached before O™ (this happens 
if 3a — 1 > (u>\/u) 2 ) 3 ), which corresponds to the case of figure and to our experimental 
parameters. The target gas can condense only if T>™ ax > 2.612, and in the same way, the 
buffer gas can condense only if D™ ax > 2.612. If these two conditions are satisfied, the 
buffer condenses first if T)^ > 2.612 otherwise the target condenses first. These three 
conditions can be written as conditions on the number N 2 of target atoms. This leads to 
define three critical numbers iV"^, iVg and iV| of target atoms, for which T)( = \ D™ ax and 
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-pmax are reS p ec tively equal to 2.612. Although each of these numbers depend on all initial 
conditions, the ratios N%/N% and N\/N% are independent of the initial values N[ m and T im 
because P™ ax /X>™ ax and X>( = ) /D™ ax do not depend on the initial conditions (see eq. |3| and 



sympathetic cooling. From the experiment, we measure N$ = 3 x 10 4 , N\ = 8 x 10 4 and 
= 2 x 10 5 which is in good agreement with the calculated values for a typical value of 
the evaporation parameter r\ ~ 6.5 (see fig. |H). 

In a second experiment, we have studied the role of thermalization in sympathetic cooling 
by changing the bias field from 56 G to 207 G. This reduces the vertical oscillation frequencies 
^Fz = {(— l) F rriF Hb,{G 2 / Bq — C) /2M} 1 / 2 where M is the mass of the atom and fiB the Bohr 
magneton, and therefore increases the relative gravitational sag A = g/ujf z — g/oj% z between 
the two atomic clouds |J . Figure |l|.b represent successive absorption images while decreasing 
the final frequency of the evaporation ramp, at a bias field of 207 G (A = 26 /mi). We initially 
observe sympathetic cooling of the target gas until a temperature of about 400 nK is reached. 
Then sympathetic cooling stops. This is in contrast with the first experiment at a bias field 
of 56 G, where the relative sag is smaller (A = 7 /xm), and sympathetic cooling works all the 
way down to BEC of the two species (see fig. [I|.a). 

To render a quantitative account of these observations, we study how the thermal contact 
between the two clouds evolves during the cooling. For that, we calculate the energy ex- 
change rate W via interspecies elastic collisions. We assume Maxwell-Boltzman distribution 
at temperatures T\ for the buffer and T 2 for the target gas (we neglect the cut-off due to 
possible evaporation). Introducing the notations vq for the center of mass velocity, v and 
v' for the relative velocities respectively before and after the collision, the energy received 
in the laboratory frame by a buffer atom during the collision is 4^ vq ■ (v — v') . At low 
temperature we only consider s-wave scattering and the term vg ■ v' averages out to zero. 
After integration over positions and velocities, we get W = fcs(^2 — T\)T, where T is the 
number of interspecies collisions per unit of time, found equal to 



|J) and because (see eq. 0). These ratios are represented in a diagram 

of the energy cutoff r\ = 3a + 2 versus the number of target atoms A^/iVJ on figure |3|. 



We can identify four regions on this diagram [T2[], corresponding to the different regimes of 
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In equation |5j, <7i 2 is the interspecies elastic cross-section (taken independent of the tempera- 



ture), V = yJk B (T! + T 2 ) /M is the RMS sum of thermal velocities, and p z = Jfyffi + 



is the RMS sum of the vertical sizes of the clouds (with similar expressions for p x and p y ). 
The exponential term in (|5|) describes the overlap of the two clouds. Thus the energy ex- 
change rate W will become vanishingly small if p z becomes smaller than the gravitational 
sag A. We can now analyze the role of the thermal contact in the experiments. At a bias 
field of 207 G (sag A = 26 /mi), and high temperature (above 400 nK), the effect of the sag 
is negligible and the clouds are thermalized. During the cooling, the sizes of the clouds and 
therefore p z decreases. At the temperature of T% ~ T 2 ~ 400 nK where we observed that 
sympathetic cooling stops, we have p z = 12 pm and the exponential term in equation |5] takes 
the value 0.1, so that the energy exchange rate is reduced by an order of magnitude. In 
contrast, at a bias field of 56 G, where A = 7 pm, the parameter p z at the temperature of 
condensation 300 nK, (p z ~ 8 pm) is still bigger than A, the collision rate T is only reduced 
by 30 % because of the overlap, and sympathetic cooling works until reaching BEG 

Starting from the interspecies energy exchange rate W, and using the total energy conser- 
vation, we can derive evolution equations for the temperatures T\ and T 2 of each species taken 
separately at thermal equilibrium. Denoting AT = T\ — T 2 , and T = (7\ + T 2 )/2, we obtain 
an interspecies thermalization rate, which, for equal trapping frequencies {uii = uj 2 = u), 
takes the simple form 

I- 1 dAT _ (Ni + N 2 ) u; 3 cr 12 M 

r ~ AT dt ~ 3k B T 2 7T 2 ' U 

This result is remarkable, by analogy with the individual thermalization rate of a single 
species r _1 ~ 7/3 []13| , where 7 = Nuj 3 aM/2ir 2 kBT is the average elastic collision rate 
in an harmonic trap. Assuming identical values for the elastic cross sections a±, a 2 , and 
o~i2 fT^! we then have r _1 ~ r-f + tvT 1 where rf/ 1 and r 2 _1 are the thermalization rates of 



each species considered separately. This means that thermalization between the two species 
happens faster than thermalization of each species. So, if the evaporation ramp is adapted 
to cooling of the buffer alone, interspecies thermalization will be efficient. 

If the masses are different, equation^ remains valid at lower order in AT/T, if we replace 
M by M. = ( 8 j(^+M2) 3 ; which is symmetric in the masses Mi of the buffer and M 2 of the target. 
The equivalent mass M. is maximum for Mi = M 2 as in our experiment, and decreases for 
uneven masses (for significant different masses, M. is much smaller than the smallest mass). 
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It is interesting to note that if equation |6| is applied to any arbitrary partition in a single 
species, where Aq and N2 = N — N± atoms are taken at different temperatures, then the 
thermalization rate between the two subsamples is exactly equal to 7/3, in agreement with 
the numerical simulations of [I3|. To our knowledge, such an analytical result had not been 
published. 

To conclude, we have first presented a simple analysis allowing one to quantitatively 
predict the outcome of sympathetic cooling, in the case of no losses in the target. We have 
then obtained an expression of the thermalization rate between the buffer and the target, 
that quantitatively describes the effect of the overlap between the two clouds. Moreover, 
when the overlap is good, the interspecies thermalization rate is the same of the intraspecies 



thermalization rate, and sympathetic cooling is very efficient || [15]. In contrast, we predict 
a strongly diminished interspecies rethermalization rate for very different masses of the 
buffer and target. Our approach allows to derive an analytical expression of the intraspecies 
thermalization rate in surprisingly good agreement with previous numerical calculation. 
Although our modeling is very simple and does not take into account quantum statistics 



effects near condensation [16], it is a useful tool, that can easily be generalized, to plan or 
analyze experiments using sympathetic cooling. 
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FIG. 1: Sympathetic cooling of 87 Rb atoms in \F = 2, mj? = 2) (upper images) by evaporatively 
cooled atoms in \F = l,mj? = —1) (lower images): effect of the thermal contact, a) Small sag 
(7/(/m) between the two clouds: sympathetic cooling yields BEC for both species, b) Larger sag 
(26 fim): sympathetic cooling stops after the third snapshot (400 nK) when decreasing final rf 
frequency. 

FIG. 2: Phase space densities of the buffer (solid line) and the target gas (dashed) as a function of 
the evaporation progress, characterized by the ratio Ni/Nf 11 of remaining buffer atoms. This plot 
corresponds to the case 3a — 1 > (wi/u^) 3 - The horizontal line is the critical value 2.612 above 
which BEC occurs. Here, the target and then the buffer will condense. 
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FIG. 3: Outcome of sympathetic cooling as a function of the number of target atoms iV^/iV^ - 
The evaporation parameter r] is the ratio of the energy cutoff to the thermal energy ksT. Trap 
parameters are fixed to u^/^i = The lines represents the critical number ratios N2/N2, N^/Nfi 
and Nrt/Nrj. They separate regions where either dual, single or no BEC can be formed. Black dots 
indicate the critical number found in our experiment. 
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